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The rupture of a medium under stress typifies breakdown phenomena. More generally, the latter 
encompass the dynamics of systems of many interacting elements governed by the interplay of a 
driving force with a pinning disorder, resulting in a macroscopic transition. A simple mean-field 
formalism incorporating these features is presented and applied to systems representative of fracture 
(3J[)' phenomena, social dilemmas, and magnets out of equilibrium. The similarities and differences in 

the corresponding mathematical structures are emphasized. The solutions are best obtained from 
a graphical method, from which very general conclusions may be drawn. In particular, the various 
classes of disorder distribution are treated without reference to a particular analytical or numerical 
form, and are found to lead to qualitatively different transitions. Finally, the notion of effective (or 
phenomenological) theory is introduced and illustrated for non- equilibrium disordered magnets. 
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When boarding passengers, airline personnel usually call a limited number of rows at a time — "Now seating rows 
20 to 30!" — to ensure an efficient and orderly filling of the aircraft. Nevertheless, there are always some hurried 
individuals who unduly join the forming line before their row is called, only to be rebuked by the flight attendant. 
In most instances the line is barely disrupted, while in some cases, such as delayed flights, the "avalanching" of 
disobedient passengers creates a long, disorderly line, absorbing nearly everyone. Why is this? Let us consider the 
mechanism by which the line forms in greater detail. First, the forty or so passengers whose rows have been called 
move toward the gate, initiating a line. Next, a few "impatients" (or "disobedients" ) add to the ranks, and perhaps 
', also some absent-minded physicists who stand up as soon as they hear an announcement and innocently approach 
£h ' the gate. Seeing this, some of the people still sitting and waiting calmly, might lose their patience and decide to join 
the growing line. This "domino effect" may die out very soon, with just a few outsiders hiding in the queue, or it 
may continue, resulting in a catastrophic boarding process. 

The line formation at a boarding gate incorporates all the elements of a breakdown phenomenon [Q : it involves a 
-H , large number of elements (individuals), each in one of several possible states (sitting or in line), which interact (the 
longer the line, the stronger the temptation), and are driven as a whole (calls of the flight attendant). Furthermore, 
different people have different "obedience thresholds" that depend on their personality, mood, schedule, etc. For a 
given length of the line, a fraction of these thresholds are crossed and the resulting disobedients stand up and move 
toward the gate. Equivalently, a sitting individual decides to stand up and join the queue with some probability 
(dependent on the line's length), thereby introducing an element of stochasticity, or disorder. Similar features are 
present in a great variety of processes, occurring, for example, in materials under tension conductors with 

an imposed current Jl],||, magnets driven by an applied field p], and social phenomena such as the "boarding-gate 
problem." The rupture of the material, destruction of the conductor, macroscopic flipping of the magnetization, or 
the failure of the social process, respectively, can be viewed as a breakdown of the system. 

The present paper aims to draw a simple quantitative picture of breakdown phenomena in disordered systems, 
unifying the various realizations into a single formalism while still allowing for the details that distinguish each case. 
The most drastic simplification underlying our approach is its "mean-field" nature, which relies on average quantities 
by ignoring fluctuations and, consequently, any idea of locality or modulated spatial correlation. Although such a 
crude description can depart significantly from three-dimensional reality, it is often accurate in higher dimensions 
and establishes a coarse framework which may then be refined to include fluctuations (perturbatively) . Nevertheless, 
. , some systems, such as magnets with long-range interactions or the boarding-gate line, are properly described by a 
mean-field theory; in the boarding gate problem, for example, any notion of locality is absent from the interactions, 
JH \ as any passenger can see the line as well as any other. 

In addition to differences intrinsic to various examples, our approach should allow for the introduction of more than 
one type of disorder, as even a given system's behavior might very well depend on the character of the stochasticity 
present. Better, we would like to divide all the possible disorder distributions into classes, each associated with a 
distinct qualitative behavior. To achieve this, we use a graphical scheme where we can draw different disorder 
distributions, enabling us to handle the latter in a very general fashion, without reference to a particular analytical or 
numerical form, and to distinguish the features relevant to the breakdown phenomenon pj. These features define the 
various "universality classes," identified by one of their (drawn) members or by a generic statement (for example, "the 
class of all continuous distributions" ) . A graphical method also avoids the possible misconceptions of a more specific 
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treatment J/J. Indeed, if the qualitative behavior of a physical system changes upon making a particular, for example, 
Gaussian, distribution wider, what is the associated key feature? Is it the variance of the Gaussian or its maximum 
value? Finally, a graphical scheme allows us to "follow" the state of the system as it is driven, thus capturing the 
whole breakdown phenomenon in a single picture. 

We develop the method in the context of fractures (Sec. II), for which it was originally devised |^^|, before applying 
it to two other realizations interesting in their own right. In Sec. Ill, some quantitative results are obtained for the 
boarding gate problem. Section IV focuses on driven magnets, a workhorse of non-equilibrium statistical mechanics 
as well as a central model for hysteresis. In Sec. V, we introduce the notion of effective (or phenomenological) theory, 
a broader framework than the models examined. We return to the discussion of driven magnets in this more general 
context by constructing a simple effective theory suited to the problem, and mention some natural extensions to 
richer, non-mcan-ficld approaches. Finally, we conclude (Sec. VI) with a summarizing picture and a brief allusion to 
dynamics. 

II. FRACTURES AND EARTHQUAKES: THE DEMOCRATIC FIBER BUNDLE MODEL 

The democratic fiber bundle model was introduced by Peirce Q in 1926, and formalized and studied in greater detail 
by Daniels || in 1945. It consists of Nq (— > oo) fibers pulled by a force F which they share uniformly (democratically). 
As the force is increased, a fiber breaks if the stress it undergoes exceeds a threshold (or strength), drawn from a 
distribution p(x). After each set of failures, the total force is redistributed over the remaining fibers. As members of 
the British Cotton Industry Research Association and the Wool Industries Research Association, respectively, Peirce 
and Daniels sought a more fundamental understanding of the 'hank' and 'lea' tests for wool and cotton yarns, in 
which a hank (or bundle) is stretched between two hooks until it ruptures. A more life-threatening illustration of the 
model is offered by the failure of an elevator rope (Fig. 1). 




FIG. 1. The democratic fiber bundle model realized by an elevator rope. 



The democratic fiber bundle model has further been applied to the study of cracks and fractures jTj and earthquakes 
H, as a simplified model that retains some of the essential features of the original problem. In a failing medium, stress 
is continually redistributed due to the nucleation, growth, and coalescence of cracks. Realistic descriptions of fractures 
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and earthquakes may lie anywhere between this democratic or infinite-range limit embodied by the democratic fiber 
bundle model, also sometimes referred to as the "global-load-sharing" or "equal-load-sharing" model, and the local or 
short-range limit, the "local-load-sharing" model H, in which the stress released by a broken bond is passed on to its 
immediate (intact) neighbors. Far from being one-dimensional or democratic, real-world systems may be (effectively) 
two- or three-dimensional, with connectivity properties that lead to a more complicated and richer picture. Alluding 
to the applicability of the model to the testing of a cloth sample rather than a hank, Daniels himself recognized || 
that "in closely woven fabrics the cross-threads afford a measure of support which introduces complications." 

Beyond the obvious question of the bundle's strength (maximum sustainable stress without full rupture) addressed 
by Peirce j2| and Daniels ||, we would like to uncover, in the large Nq limit, its behavior characterized by the number 
N(F) of intact fibers as the force F increases from to oo. In particular, we investigate the critical properties of the 
fiber bundle upon approaching a macroscopic failure, and the generic classes in which they fall. 

Upon incrementing the force from F — dF to F, a few fibers whose thresholds have been overshot, break. As a 
result, the stress increases on the remaining intact fibers, leading to a secondary set of failures. This triggers a tertiary 
set of failures, and so on, until the bundle stabilizes at a non- vanishing N(F) or snaps off as a whole (N(F) = 0). 
Consider the zth rupture which leaves Ni unbroken fibers, each now under a stress F/Ni. The thresholds of a fraction 

of intact fibers have been exceeded, namely Jp/j^ 1 p(x)dx, and the next set of failures brings their number down to 



F/Ni 



* F/Ni. 

where P is the cumulative distribution defined by 



N l+1 =Ni-N p(x)dx =N il-P[ — )}, (1) 



F 



N, 



P(x) = / p(x')dx'. (2) 



i) 



Clearly, the function N(F), defined as the number of fibers remaining unbroken under a total force F, is none other 
than Noo calculated at F (the limiting value of Ni as i — + oo). 

Problems stated in terms of a recursion relation often have an elegant graphical solution which displays in a single 
picture the presence or absence of convergence. Here, a graphical scheme for the iteration is best constructed with 
the quantities ir, = N/F, f — F/Nq, and 

*(x) = l-p(±y (3) 

in terms of which Eq. (|l|) becomes 

fxi+i = n(xi). (4) 

The function tt(x) increases monotonically from at x = to 1 at x — oo and a graphical iteration of Eq. (|]) shows 
that N(F) = Nso is given by the right-most intersection of the curve y = ir(x) with the straight line y = fx. As the 
force is increased, the straight line becomes steeper and the abscissa N(F)/F of the intersection point moves to the 
left. Figure 2 illustrates the graphical scheme for three successive increments of the force. First, / is switched from 
to fa for /—►(), the straight line is flat and intersects tt(x) at x — oo (as expected from xo = N /F) and y = 1. After 
the sudden switch, fibers with a threshold smaller than fx immediately break, causing more ruptures according to 
Eq. (^J). The corresponding iterations are represented by dotted arrows, and terminate at the intersection point with 
abscissa Xoo (fx)- Next, / is switched from fx to /2, and similarly, Xi converges to x^ (f%)- Finally, / is switched to 
/3, and the iterations lead all the way to x = 0, which corresponds to a fully ruptured bundle. If the force, rather than 
suddenly switching from fx to /2, goes through several intermediary steps, x successively converges to the abscissa 
of the corresponding intersection points, but again ends up at x (X) (f2) when / = fi. Thus N(F) is a well defined 
function of F, independent of the past history of the force as long as it is increased monotonically. 



3 







X 



x-(fi) 



FIG. 2. Illustration of the graphical scheme for three successive increments in the force. The dotted arrows represent 
iterations of Eq. (^). The latter terminate, for each value of the force, at the right-most intersection, yielding the number of 
intact fibers as N(F) = Fxca- 



In what follows, we focus on a continuous steady increase of F, which can be thought of as the limit of many small 
successive increments. We assume that the iterations occur fast enough, or equivalently, that F varies slowly enough 
for the convergence to be considered immediate. By analogy to thcrmodynamical processes, we may say that the 
force is increased quasi- statically. 

We are now in a position to investigate the critical behavior of the democratic fiber bundle model with various 
threshold distributions p(x). The graphical method provides generic conclusions, based on a (drawn) arbitrary curve 
y = 7r(x). Thus the results are independent of a specific analytical or numerical form for p(x), but follow from the 
qualitative features of the curve. The latter identify a class of distributions, associated with a given behavior of the 
bundle. For continuous distributions, there are three generic classes as illustrated in Fig. 3. These classes are 
distinguished by the character of p (x) at large argument, which dictates the bundle's behavior under large forces, 
through the slope 



d,Tr(x) 




(5) 
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FIG. 3. The democratic fiber bundle model with continuous threshold distributions, (a) Illustration of the graphical scheme. 
As the force is increased, the straight line becomes steeper and the abscissa of the intersection points, proportional to the 
number of intact fibers, move to the left. Typical solutions are shown in (b). Intact fibers are always left in bundles of type (i), 
which never rupture completely. Bundles of type (ii) fail continuously at a force i*c(«) = ^oQ, whereas bundles of type (Hi) 
fail discontinuously at F C (ijj). 



Class (i). For p(x) ~ x~~ r (x — > oo) with r < 2 (normalization requires r > 1), that is, distributions with an 
appreciable fraction of robust fibers, the bundle never fails completely. Under any force F, some intact fibers are left. 
Class (ii). In the limiting case where p(x) = ax~ 2 as x — > oo, the number of intact fibers goes continuously to zero 
at F c = Noa. Furthermore, the response function of the bundle or breaking rate dN/dF diverges as (F c — F) 1 ^ s , 
where s is the subleading power entering the distribution function, p (x) = ax~ 2 + [3x~ s + . . . , s > 2. Following 
the nomenclature adopted in the field of critical phenomena, we may refer to the macroscopic rupture of the bundle 
associated with a diverging response function as a second-order transition with exponent 7 = 1/s. 

An example of a "fat-tail distribution" is p(x) — (a/x r ) e~ a / xT 1 for any a, with 1 < r < 2 in class (i) and r = 2 
(and exponent 7 = (2r — 1) = 1/3) in class (ii). 

Class (Hi). For any narrower distribution, such that x 2 p(x) — > for x — ► 00, there is a macroscopic failure of the 
bundle at a finite critical force F c at which the number of intact fibers falls discontinuously from N(F C ) to 0. The 
breaking rate diverges in general as 

A "trivial" exponent of 7 = —1/2 is to be expected within a mean-field approach which ignores the effect of spatial 
fluctuations (see Sec. V). Of course, it is possible to choose a particular distribution with a "bump" such that 7 = l/2n, 
with n an integer, instead of 7 = 1/2, but it reflects a fine-tuned choice rather than the generic one. All continuous 
distributions with finite mean, subsumed in class (Hi), lead to a macroscopic failure of the bundle at a finite critical 
force, signaled by a diverging breaking rate. Although this divergence in the response is still reminiscent of a second- 
order transition, the situation differs from the usual case where the order parameter (here N(F)) is continuous at the 
singularity. To avoid any confusion, we refer to ruptures of type (Hi) as soft failures. 

Two examples of distributions which belong to class (Hi) are p(x) = (x/\ 2 ) e~ x / x and p (x) — (x/X) e~ x l 2X . Other 
examples include finite-support distributions, which are non- vanishing only in an interval X — \xq,xq + A], provided 
p(x) is continuous on X and p(xq) = p(xq + X) = 0. 

Fig. 3(b) illustrates the generic histories N(F) for each of the three classes. Clearly, by having more "bumps" 
(n( m \x) on Fig. 3(a) has one), that is, variations in the convexity of n(x), the bundle may undergo additional 
macroscopic (finite fraction of A^o) failures at smaller forces. None of these, however, may lead to a complete rupture, 
which always corresponds to one of the three classes described above. 
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FIG. 4. The democratic fiber bundle model for a markedly discontinuous distribution of the thresholds. In (a), the graphical 
scheme leads to (b), an abrupt failure of the bundle. 

The examples given for class (Hi) are valid for any value of the parameter A. Thus, the threshold distribution 
may be as narrow and peaked (A — * 0) as desired and still lead to a soft failure, provided p(x) is continuous. An 
abrupt rupture of the bundle with no divergence in the breaking rate preceding it is possible only if ir(x) has a non- 
diffcrentiable point (Fig. 4) Q], which in turn requires a discontinuity in p. Thus, abrupt failures in this mean-field 
model, far from being a generic feature associated with narrow disorder distributions, are an artifact of a singular 
point in p. The phenomenon is in fact even more particular, because the magnitude of the discontinuity matters. In 
the case of a finite-support distribution on X with p(xo) 7^ 0, for example, no fiber breaks up to a force Fq = NqXq. 
At F , an abrupt failure occurs only if dir (1/xq) jdx > Xq, as can be shown by the graphical method, which requires 
a minimal jump p(xq) > 1/xq. Whether this macroscopic failure yields full rupture of the bundle or not depends on 
the detailed form of p. If it does not, and if p is continuous on X, the global bundle failure is soft and belongs to 
one of the three cases examined above (case (Hi) if (x) = J dx x p(x) < oo). The description of a large Nq bundle 
by a continuous distribution p would seem more physical, unless it is made of fibers of a few different types, each 
with a well denned threshold. Consider for example a "bimodal" bundle with pN$ (0 < p < 1) fibers of strength /i 
and (1 — p) Nq fibers of strength fa < fi ||. It is easy to see from the graphical method (Fig. 5) that if /2//1 > p, 
the entire bundle snaps off abruptly at Fi — NqJi- On the other hand, if /2//1 < p it undergoes an abrupt failure 
at F2 leaving the pNo stronger fibers intact, and then ruptures completely (again abruptly) when the force reaches 
Fi = N fv 
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FIG. 5. The democratic fiber bundle model for a "bimodal" distribution consisting in pNo fibers of strength /i and (I — p) iVo 
fibers of strength /2. (a) Illustration of the graphical scheme and (b) typical solutions. The bundle fails in either one or two 
abrupt steps. 



A physical medium which fails as a result of the growth of a single large crack is often referred to as "ductile," in 
analogy with flexible materials, such as most metals, which rearrange under stress and can thus be pulled into long 
thin threads. At the other end of the spectrum, fragile or "brittle" media are characterized by a sudden failure due 
to the coalescence of many microcracks. In the context of our mean-field model, we may identify abrupt failures from 
markedly discontinuous distributions (such as narrow finite-support ones or bimodal ones) with brittleness, and soft 
failures from the continuous distributions of classes (ii) and (m) with ductility. Thus the character of a medium's 
strength depends on the shape of the distribution of randomness. 

The failure of an element, whether a fiber in an elevator rope or a patch of a tectonic plate, is often accompanied by 
the emission of a sound [|10| . Before an abrupt failure, nothing unusual is heard, whereas an increasingly loud noise 
is produced by a soft failure. Thus impurities, dirt, and disorder might be salutary in times of calamity! 



—Mi < P 

/2//l>P 
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III. SOCIAL LOGISTICS: THE BOARDING GATE PROBLEM 



Armed with the method developed in Sec II, we return to the boarding-gate problem of the introduction, and treat 
it in a quantitative fashion. The similarity between the breaking of a fiber bundle and the growth of a boarding-gate 
line becomes more apparent through the following mapping of the degrees of freedom, disorder, driving agent, and 
coupling: 

sitting person, waiting < — ► intact fiber, 

standing person, in line < — ► broken fiber, 

"obedience threshold" < — ► strength threshold, 
flight attendant's calls < — ► driving force, 

temptation from seeing the line < — > redistribution of the force. 

Nevertheless, the two problems are different in that the stress applied to any one intact fiber is inversely proportional 
to their total number, whereas the temptation here grows with the length of the line. Furthermore, for simplicity we 
focus on a single call of the flight attendant, say the first one. In the language of the democratic fiber bundle model, 
this corresponds to switching the force from zero to some value, and leaving it there. Finally, the "thermodynamic 
limit" we treat (in implicitly assuming infinitely many passengers) may be a poor approximation for a situation 
involving a few dozen or a few hundred people. The fact that the first group of passengers called, and consequently 
the ones still waiting, may not be representative of the overall threshold distribution, for example, is one of the 
possibly important finite-size effects neglected. 

Let us say that the flight attendant first calls £ passengers, who readily form a line. Seeing the latter, a few 
disobedients join in, equal in number to (JV — £)P(£), where iV is the total number of passengers waiting to board 
the flight and P is the cumulative threshold distribution. More precisely, p(x)dx = (dP/dx)dx is the probability that 
an individual unduly joins the line, when it reaches a length between x and x + dx. The process might stop at this 
point, leaving the few intruders in line, or might keep on "avalanching." In general, we can write a recursive relation 
describing the growth of the line, analogous to Eq. (|l|), as 

n i+1 =£+(N-£)P(ni), (7) 

where rij is the number of people in line after i iterations (no = £)■ The avalanching terminates at the left-most 
intersection of the line (n — £)/ (N — £) with the curve P(n) and the eventual length of the line, rioo , is read off as the 
abscissa of the intersection point. In other words, the length of the line is obtained as the smallest solution to 

" ~ ' P(n), (8) 



N -I 



leading to the following possible scenarios. 

• Most optimistically, no one's threshold is lower than the length of the initial line, that is, p(x < I) = (Fig. 6). No 
unexpected passenger shows up in the line, and the boarding runs smoothly. 

• More generally, for a large airplane with hundreds of passengers, we expect thresholds to be widely distributed. 
We consider distributions with p(x) ^ for all x £ [0,N] and P(N) — 1, according to which there arc always some 
disobedients who join the ranks no matter how small the temptation and no one is restrained enough to stay seated 
when the line includes nearly all the passengers. In this case, the behavior of the line's growth is governed by the 
character of p at large argument (Fig. 7), reminiscent of the solution of the democratic fiber bundle model with 
continuous distributions. If p(N) < l/N, the line may anarchically absorb all passengers. If p(N) > 1/N, on the 
other hand, it must stabilize at £ < < N, the precise value of which depends on the detailed form of the function 
P- 
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FIG. 6. Graphical scheme for the boarding-gate problem with a threshold distribution vanishing below I, The solid and 
dashed curves represent two possible cases, both leading to lines with no intruders. 




FIG. 7. Graphical scheme for the boarding-gate problem with non- vanishing distributions. If p(N) > 1/N (solid line), P(n) 
must cross the straight line (n — £) / (N — £) at some n smaller than N. If p(N) < 1/N (dashed lines), it may or may not. 



• Often, there is a small fraction n p /N of complete disobedients, a small fraction n s /N of complete obedients, and 
the rest of the passengers' thresholds scattered in between. For such a contingency, the distribution reads 

p(x) = ^S(x)+g(x) + ^r6(x-N), (9) 

where g(x) is a smooth function which, in a realistic model, is appreciable between two bounds m and n^. Then, the 
growth of the line saturates at a length £ + (-^rO n p < < N — f^^) n s, which depends on the relative values 
of £, n p , and n\ (Fig. 8). If n\ > £ + n p , the eventual line includes few unexpected passengers except the physicists 
— a situation easy to cope with, due to the small number of physicists. However, if n\ < £ + n p (and in particular 
if n,2 — ri\ is small (g very peaked)), the few physicists drag flocks of others along, resulting in a line with all the 
passengers but the complete obedients and perhaps a few highly obedient individuals. 
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FIG. 8. Graphical scheme for the boarding-gate problem with a distribution reflecting a fraction of complete disobedients, 
a fraction of complete obedients, and the rest of the passengers peaked in between. The character of the solution depends on 
whether the peak is closer to the disobedient (solid line) or to the obedient (dashed line) side. 



There is no doubt that social phenomena involve a great many parameters making them highly complex in compar- 
ison to "simple" physical systems such as a fiber bundle. However, in much the same way as the breaking of a wool 
hank is faithfully described without recourse to the molecular structure of the individual threads, a coarse picture of 
some human phenomena may be drawn while ignoring intricate psychological and other processes that each individual 
goes through. And it is amusing that such a phenomenological approach yields a number of quantitative, testable 
results for the behavior of a human group. 

IV. MAGNETS OUT OF EQUILIBRIUM: CRITICAL HYSTERESIS IN A DRIVEN ISING MODEL 

In the most modest speck of dust, an incredible number of electrons are dancing around, creating small currents 
and with them magnetic moments. Each electron also carries its own spontaneous moment, and it is the complicated 
addition of all these effects which determines the total magnetization of the dust. In the spirit of foregoing intricate 
details in favor of a coarser but simpler picture, Ising and Lenz proposed a model |llf| in which a material is divided 
into many elements (labeled by their position r), each carrying a (unit) magnetic moment, or spin, which tends to 
align with the local magnetic field it experiences, like a little compass. In the so-called "Ising model," the spins 
are constrained to a given direction, may point upward (m(r) = +1) or downward (m(r) = +1). The average 
magnetization is defined as 



where N — > oo is the total number of spins, and the local field is composed of the applied field and the one created by 
neighboring spins. On a mean-field level, one assumes that the system may still be faithfully described if the many 
degrees of freedom are reduced to a single quantity, here the average magnetization. This assumption becomes exact 
if the interaction between spins is infinite in range, that is, all spins interact with all others in the same way, so that 
they contribute to the local field equally everywhere. Even if the coupling has a limited range, we can imagine that 
such an approach becomes more and more accurate as the dimensionality is increased, as spins acquire more and more 
neighbors, thus "averaging out" their contributions to the local field. In fact, the mean-field approximation yields the 
correct critical behavior above a given, finite dimension (see Sec. V). 
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FIG. 9. Illustration of a few "sites" of a two-dimensional Ising model with a random field. The thick arrows represent the 
spins which can point upward (m = +1) or downward (m = —1), and the thin arrows represent the random field h. Each spin 
experiences a uniform applied field H . 



We confine our study to zero temperature, so that spins are always perfectly aligned with the local field, forbidding 
any thermal fluctuation. In real materials, however, the presence of impurities leads to random fluctuations. In many 
cases, the effect of the impurities is to perturb the local magnetic field which acquires a stochastic component in 
addition to a uniform one (Fig. 9) rfTij . More precisely, the mean- field, driven, disordered version of the Ising model, 
first introduced in Ref. ^, is defined by a local field 

/(r) = H + h(r) + JM, (11) 

where H is a uniform field (chosen by the experimentalist), h is a random field chosen at each position r from some 
probability distribution p(h), and J > represents the strength of the spin-spin coupling. The configuration of the 
system is then given by the "alignment rule" 

m(r) = -1 if /(r) < & h(r) < -H - JM, . , 

m(r) = +1 if /(r) > & h(r) > -H - JM. ^ ' 

How does the driven disordered Ising model evolve when the (driving) field H goes from — oo to oo? Initially 
(H = — oo), all the spins point downward (m = —1). As H is increased from — oo, the local fields with large positive 
random fields soon become positive, and the associated spins flip upward. These flips further increases the local field, 
which induces new spin flips, and so on. Once H is incremented enough, the negative h's tend to pin the corresponding 
spins downward, and the configuration results from the competition between the driving agent H, the "restoring" 
force JM, and the pinning forces h — a common scenario for forced or driven random media p5| , ^6| . 

The avalanching present in this simple model is reminiscent of fracture phenomena and the driven disordered 
Ising model can be mapped to the democratic fiber bundle model as 

downward spin (m = —1) < — ► intact fiber, 

upward spin (m = +1) < — ► broken fiber, 

random field (h) < — ► strength threshold, 

uniform field (H) < — > driving force, 

spin-spin coupling (J) < — > redistribution of the force. 

As in the case of the boarding gate problem, the main difference with the democratic fiber bundle model arises from 
the way in which the individual degrees of freedom interact, and the iterative equation for the fraction of downward 
spins, in the driven disordered Ising model, reads 

n ii+1 =P(-H- JMi), (13) 

where P(x) = p(h)dh is the cumulative distribution function. Equivalently, in terms of the magnetization 
M = l-2n ; , 
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M l+1 = 1 - 2P(-H - JMi). 



(14) 



As before, the right-hand-side is a monotonically increasing function of the avalanching variable M^, and the so-called 
hysteresis curve M = M QO (H) is given by the left-most intersection of y(M) = 1 — 2P(—H — JM) with the bisector 
(y = x). 

Rather than considering all possible types of disorders, we focus on three, commonly encountered in the literature 

(1) Gaussian disorder. p(h(r)) — -j===er h l 2a . As H is increased, the curve y{M) is translated to the left and 
the intersection point consequently moves up (and to the right). Three qualitatively different solutions are obtained, 
depending on whether the maximum slope y'(M = —H/J) = 2 J j\j2-Ka 2 is greater, equal, or lesser than 1 (Fig. 10). 





FIG. 10. The driven Ising model with a Gaussian random field, (a) Illustration of the graphical scheme for sub-critical, 
critical, and super-critical values of the disorder. The curves y(M) are drawn at H — 0; as H is increased, they move to the 
left, starting in the far right at large negative H and ending in the far left at large positive H. (b) Corresponding typical 
magnetization, or hysteresis, curves. The dashed line represents the case of a decreasing H , in the sub-critical regime. In the 
critical and super-critical regimes, the curves for increasing and decreasing H are superposed. 
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• fj < a c = For a small disorder, there is a jump in the magnetization at some field -ffjump > 0), signaled, as 
in the democratic fiber bundle model, by a square root divergence of the response function or susceptibility 

X=^-(H iump -H)-^. (15) 

Larger values of a reduce the magnitude of the discontinuity, which ultimately shrinks to a point at a — \j2jn3 . 

• a = a c : A specific value a c of the disorder yields a critical hysteresis curve, which, although continuous, displays a 
singular point at H = M = 0. Expanding y(M), we easily obtain the nature of the singularity, as 

X(H -> 0) - H- 2/3 . (16) 

In a more general line of reasoning (see Sec. V), we note that the asymmetry of y about the origin (at H = 0) rules 
out even powers relating M to H as 

H ~ M 2n + higher powers, (17) 

and in particular a quadratic dependence H ~ M , leaving a cubic behavior H ~ M 3 as the best candidate; inverting 
and differentiating readily leads to Eq. (|l6|). 

• a > <r c : At large disorders, the magnetization curve is smooth, with a finite susceptibility everywhere. 

Figure 10(b) shows generic curves corresponding to the above three cases. This phase diagram is very similar to 
that of the democratic fiber bundle model (classes (i)— (in)), in that a discontinuous, first-order-like region terminates 
at a critical, or second-order, point with a different (non-trivial) exponent. In fact, it is believed [n9U2Q] that, while 
the global phase diagram is unchanged in a local theory extending beyond a mean-field approach, the precursor 
divergence of Eq. ( jig ) is suppressed, resulting in proper first-order behavior. (If a similar conclusion applies to 
fracture phenomena, it might help to explain the presumably abrupt failures often observed.) Furthermore, as a 
approaches cr c from below, -f/jump decreases to zero, progressively hampering the hysteretic (lagging) nature of the 
solution; for a > er c , both the remnant and coercive fields are vanishing. Again, this is an artifact [M of the present 
formulation. These remarks point at the tip of an iceberg onto which the mean-field approximation often runs. 
(2) Uniform disorder. p(h(v)) = 1/2(7 for — a < h < a and vanishes elsewhere. This form may be viewed as a 
"digitized," roughened, or sharpened version of Gaussian disorder. Again, the character of the solution depends on 
the magnitude of the slope y'(M = —H/J) = J jo relative to 1 (Fig. 11). 
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FIG. 11. The driven Ising model with a uniform random field, (a) Illustration of the graphical scheme and (b) typical 
solutions. 

• If a < J (small disorder), the magnetization switches from —1 to +1 at -ffj ump = J — cr > 0. 

• If a > J (large disorder), the magnetization is continuous and piecewise linear; it starts to increase linearly from 
— 1 at H = J — <t, and reaches +1 at H = a — J. 

As in the previous examples, the nature of the solution varies significantly from one type of disorder to the next. 
In particular, we note that the discontinuity of the uniform distribution rules out any divergence or critical point, in 
agreement with our observations for the democratic fiber bundle model. 

(3) Bimodal disorder. p{h{v)) = [8 (h + cr) + 5(h — a)] /2. The bimodal distribution is often used as an alternative 
to the Gaussian distribution, notably in computer simulations. It represents an even more "digitized," or sharpened, 
version of disorder, reminiscent of the democratic fiber bundle model with fibers of two different types. The graphical 
scheme is illustrated in Fig. 12, along with typical solutions which fall, as above, into one of two classes. The phase 
diagram is especially simple to understand for a bimodal disorder, which limits random fluctuations as much as 
possible and reduces the problem to a two-body one. Half of the spins (those with h = +cr) flip upward when H 
reaches -f/jumpi = J — cr, leading to a vanishing magnetization. This further triggers the flipping of the remaining 
spins if their (updated) local field / = -ffjumpi — o~ is positive, that is, if a < J/2 (small disorder), resulting in a 
global switch from M = — 1 to M = +1 at -ffjumpi- If cr > J/2 (large disorder), on the other hand, the magnetization 
remains at zero until the local field changes sign, at i/j U mp2 = cr, when the remaining downward spins flip to yield 



M = +1. 
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FIG. 12. The driven Ising model with a bimodal random field, (a) Illustration of the graphical scheme and (b) typical 
solutions. 



In contrast to the Gaussian and uniform cases, the large-disorder phase exhibits hysteresis: M(H) is not symmetric 
about H — 0, there is a remnant magnetization M = — 1 (if -ffjumpi > that is, a < J) and a coercive field which can 
be defined as -ffj U mp2- Clearly, a hysteretic behavior originates in the coexistence of more than one possible solution 
to a dynamical system [|2l|| . Rather than being chosen by an optimum principle, the actual, physical solution follows 
from the history, or, in the case of a monotonically increasing H, from the initial conditions. Thus, metastability 
and hysteresis are intimately related. The random fluctuations tend to pin the system to a metastable state. The 
thermal fluctuations, which we have ignored, play an opposite part in sweeping the system in its configuration space, 
thus favoring ergodicity. An important question relates to the interplay of these two effects. It is believed that 
disorder generally dominates and that thermal fluctuations are irrelevant at criticality, both in equilibrium []l4f and 
in non-equilibrium [p 19 2l[] problems. 



V. DIGRESSION: EFFECTIVE THEORIES 



The path from elementary laws to the multitude of fascinating observed phenomena is definitely not a straight and 
plain one. The situation is a bit like that of a shadow show in which shadows of dragons, wolves, and birds, are cast 
on a white drape by the moving hands of a pantomime. The goal is to understand the unfolding story. Focusing on 
the artist's hands is probably not the most fruitful tactic, and surely not the easiest. The shadows allow a more direct 
grasp of the story by retaining the important aspects of the hands' intricate movements and shapes, while smearing 
away irrelevant details such as the nails or the hair on the fingers. In theoretical physics, the mime's hands are the 
fundamental laws, and the shadows are effective or phenomenological theories. 

Hydrodynamics constitutes a prominent example of an effective theory, where the molecular nature of a fluid is 
disregarded in favor of a description in terms of a homogeneous medium. The collective motion of hundreds or millions 
of particles is reduced to that of their center of mass, that is, to a local average, which smears out a large number 
of details and fine characteristics, leaving a framework with a smaller number of variables than the original one. 
Typically, a fluid is then described by fields of density, velocity vector, and temperature — five quantities associated 
with each point in space and time substitute for the many molecular degrees of freedom around that point. 

A delicate question relates to how much of the finer print, or fluctuations, needs to be included for an accurate 
description of a given phenomenon at a given length scale. Various effective theories may be constructed to describe 
a given system, the ones with the fewest allowed fluctuations being termed "mean-field." Often, spatial fluctuations 
are ignored altogether resulting in highly coherent behavior, and the system reduces to a single (or a few) degrees of 
freedom. 

The hydrodynamic equations can be derived from symmetry and conservation principles alone. These principles 
are satisfied at the molecular level and unspoiled by any averaging procedure. Thus, an effective theory ensues from 
robust properties that span many length scales, and is free from a specific (mechanistic) model. It has lost the details 
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involved in the latter, but gained the sturdiness of generality: an effective formulation encompasses all possible models 
lying within the constraints imposed by a set of symmetries. 

It is therefore worthwhile to construct a phenomenological theory of breakdown phenomena; in what follows, we 
shall concentrate on the driven disordered Ising model [ppO[|. In addition to its robustness, an effective field theory 
offers an analytical framework in which the range of the coupling between spins, the dimensionality of the systems, or 
other aspects may be modified and studied in a unified fashion, and for which a number of powerful tools have been 
developed f22fl . 

In the spirit of hydrodynamics, we consider a spin field or density s(r) regarded as a local average of magnetic 
moments at r € 3? d , where d is the dimension of the system. Its value lies anywhere within an interval centered at 
zero, and we further generalize the problem by allowing s(r) S (—00,00). Inspired by the model [p|j20| of Sec. IV and 
its solution as given by y(M) — M = 0, we posit that the spin density at each point is obtained from an algebraic 
equation 

f r ({s(v)},H,h(r))=0 (for all r), (18) 

where H is a uniform applied magnetic field and h(r) is a Gaussian random field felt by the spins at r. This system 
of equations may have more than one solution; if H is monotonically increased as before, the spins start off very 
negative, and the physical solution corresponds to the "left-most" one, with smallest (most negative) s(r)'s. If, 
rather, H is monotonically decreased from a high positive value, we expect the same set of solutions with an opposite 
sign, consistent with the reflection symmetry 

fr({s(r)}, -H, -h(r)) = f T ({s(r)}, H, h(r)), (19) 

a first constraint on the form of / r . As before, the quantities of interest are the overall spatial average (s) of s(r) 
and its response to the applied field. An effective, coarse-grained theory is clearly best suited to answer questions 
on the system's global properties, and in particular on its critical behavior which involves a collective response to 
an infinitesimal change in H. From the study of the driven disordered Ising model, we expect the critical point to 
occur at H — (s) — 0. (We note that this is also a special point in terms of the symmetries, due to the statistical 
invariance under h — > —h.) Close to the latter, we may describe / r by an expansion in its arguments. In contrast to 
fluid dynamics or pure critical phenomena, here quenched disorder imposes some spatial fluctuations, which we limit 
by considering a small random field, leading to 



fr = 



J d d r' J(r, r')s(r') + AH + Bh(r) 

+ J d d r'd d r"d d r"'K(r, r', r", r"') s(r')s(r")s(r"') 

+ hJ d d r'd d r"L(r, r', r") s(r')s(r") + h(r) J d d r'd d r"M(r, r', r") s(r's(r") 

+ H 2 J d d r'N(r, r') s(r') + h(r) 2 J d d r'0(r, r') s(r') + Hh(r) J d d r'P(r, r') s(r') 



CH 3 + Dh(r) 



Along the lines of Eq. (|20|) , where we have written the most general expansion consistent with the governing symme- 
tries, we further expand the interactions. For example, we write 

J dVj(r,r>(r') = J dVj(r,r') | S (r) + £ (r Q - r' a ) A a ( r ) 

a, P H ) 

= || dVj(r,r')| s(v) + E {/ dV ^ - r '«) J M} ^rs(v) 
+ k E {/ dV ^ r ' a ) h ~ r 'p) J ( r > r ')| 



d 2 

s(r) + ... (20) 



drpdr. 

If the couplings are invariant under translations and rotations, then J(r,r') = J(|r — r'|), the odd powers in Eq. ( |20| ) 
vanish and the even powers take a symmetric form, resulting in 
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d d r' J(r, r')s(r') = as(r) + 6V 2 s(r) + . . . , 



(21) 



where a and b are numerical coefficients obtained from J. The derivative terms couple neighboring spins by hampering 
(or favoring, depending on the sign of the prefactors) spatial fluctuations. In their interplay with disorder effects, 
such short-range interactions introduce higher complications. By analogy with the driven disordered Ising model, 
we consider a mean-field approximation by substituting the derivative interactions with a uniform restoring force 
composed of (s), (s) , H (s) , H 2 (s), etc. If a similar procedure is applied to the various terms in / r , we obtain 

/ r = as(r) + b' (s) + AH + Bh(r) + cs(r) 3 

+ pother terms, such as (s) 3 , Hs 2 , h 2 s , or H 3 ^j . (22) 

To get a flavor of the construction and treatment of a phenomenological theory, rather than to examine the problem 
in full rigor, we abandon the "other terms" for now — we argue below that they do not change the critical properties. 
Furthermore, for ease of notation, we shall not drag along the many prefactors; the reader may either imagine them 
present or assume that they have been set to unity. Thus, our task is now to solve Eq. ([l8|), which has been reduced 
to 

s(r) + (s) - H + h(r) + s(r) 3 = 0. (23) 

The prefactor of H has been set to —1 so that s aligns along H at large values of the latter. This stability condition 
is imposed in addition to the various symmetries. In solving Eq. (p3|), we think of (s) as a given function of H. If 
the resulting s(r) does not average to (s), our choice was not good, and we start over with a better guess for (s), 
until eventually the spatial average of s(r) is equal to (s) for all H. In practice, Eq. (^3| ) may be solved iteratively by 
reshuffling the terms, as 

s(r) = - (s) + H - h{r) - s{r) 3 

= -{s)+H- h(r) - [- (s) +H- h(r)} 3 (24) 

Clearly, performing the average over space or over disorder is equivalent, and the self-consistency condition results in 

(3cr 2 - 2) (s) + (s) 3 = (3a 2 - 1)H + (higher-order terms) , (25) 

where a 2 — (h 2 '). For large disorders, the magnetization (s) is a non-singular function of H. At a critical value 
a 2 = 2/3 = er 2 , however, (s) cx H 1 / 3 and the response function diverges as 

~ £T 2/3 , (26) 

reproducing the critical exponent of the driven disordered Ising model. 

Let us return to the "other terms" neglected in Eq. (p3|), and argue that they do not modify the exponents. Indeed, 
any one of them either vanishes upon averaging (h 3 for example) or is subleading in the limit H, (s) — > (for example, 
H 2 s yields H 3 and H 2 (s) upon averaging). In a third possibility, the additional terms modify the prefactors in 
Eq. (|2^) (for example, h 2 s yields a 2 (s) upon averaging), but do not generate new types of terms. The symmetries 
are constraining enough to impose the cubic root behavior of the magnetization close to criticality. In a mean-field 
approach, in particular, where the system is reduced to a single degree of freedom, the exponents follow from an 
algebraic equation such as Eq. (p5|). This observation sheds some light on the (then quite cryptic) remark in Sec. II, 
that an exponent —1/2 "is to be expected within a mean-field approach." Without the reflection symmetry that we 
have assumed in the derivation of Eq. (|2^) , quadratic terms would appear and clearly lead to an inverse square root 
divergence in the response function. 

Another way to understand the decoupling that occurs in a mean-field theory and its essentially zero-dimensional 
nature (single degree of freedom) is obtained in a formulation closer to equilibrium statistical mechanics where the 



central objects are probability weights and partition functions. From rewriting Eq. ( p3[ ) as 

- h(r) = s(r) + (s) - H + s{r) 3 (27) 

and the random field density 

p(h) = We~2&, (28) 



dH 
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where J\f is a normalizing factor, the spin field is distributed according to 



p (s(r)) = A/"exp 



2a 2 

Again, we require the self-consistency condition 



-L ( s(r) + ( S ) _ H + s ( r )3) 



x ||Jacobian|| . (29) 



(s) = J ds s p(s), (30) 

and the problem reduces to a zero dimensional one: finding the average of s at a single "site." The right-hand-side of 
Eq. (30) may be evaluated by expanding the exponential in the super-quadratic terms of its argument and calculating 
Gaussian integrals. 

In the absence of modulated interactions, the Jacobian is easily evaluated as 1 + 3s. Beyond the mean-field picture, 
however, fluctuations matter regardless of their coupling to h, and the need for and use of a Jacobian becomes quite 
cumbersome. For a brief (and very superficial) discussion of richer theories, we return to the formulation of Eq. (pq), 
whose analog for the case of short-range interactions reads 

s(r) + V 2 s(r) - H + h{r) + s(r) 3 = 0. (31) 

How does the coupling between neighboring spins, expressed in the Laplacian term, modify the critical behavior? 
This question is best examined through the renormalization-group theory [p4f , which provides a systematic way to 
keep track of the role of each intervening term. Roughly speaking, it enables us to carry out an iterative solution (or 
averaging), similar to that of Eq. (pi|), in successive steps corresponding to progressively larger scales. Interestingly, 
it is found that the answer depends only on the dimension d of the system. Above six dimensions, a system described 



by Eq. (31) behaves just as the mean-field one of Eq. ( (23|) |19| , [20|| , as far as critical exponents go. In some sense, it 
is bundled tightly enough so that a small push might move a lot. (Many hay stems are often dragged along when 
one is pulled out of a three-dimensional pile; this does not happen if the hay stems are laying on the two dimensional 
surface of a pond.) In d < 6 = d c (called "upper critical dimension"), the situation is far less trivial, and novel 
exponents emerge from the interplay of the Laplacian and random terms ]l9| , p0[ . We might also expect intuitively 
that longer ranged interactions "more easily" lead to mean-field behavior and, indeed, they reduce the upper critical 
dimension. The case in which d c < 3 explains the trivial, mean-field-like exponents associated with some observed 
critical phenomena that, a priori, involve many complicated interactions. 

In some instances, fluctuations bring about further surprises. A mean-field approximation somehow "homogenizes" 
the system, imposing a high degree of symmetry. Fluctuations might then break the latter into lower symmetries: 
a flat surface is invariant under any translation, whereas a periodically corrugated one is unchanged only under 
translations by the period's multiples. New critical points then appear, corresponding to the new symmetries. This 



is exactly what happens in a system defined by the vectorial generalization of Eq. ( |23] ) |21|, with the new degrees of 
freedom s(r) e P. At the mean-field level, the critical properties correspond to a full rotational symmetry (in spin 
space). For short-range interactions, new critical points appear, associated with rotational invariance in subspaces of 
5J" llj- 
lt is truly remarkable that we are able to make precise predictions based on such a general approach, without a 
definite model or "fundamental" understanding of the details. It is equally surprising that a complicated assembly of 
electrons and nuclei generates pure numbers, the critical exponents, which can be measured with great accuracy. 



VI. CONCLUSION 



The models examined in Sees II-IV are probably less interesting for their own sake than for their illustrative power. 
Also, they provide motivation and inspiration for more complex, robust formulations closer to physical reality. This is 
especially true of simple models involving ubiquitous features such as driving forces, elastic interactions, and disorder. 
The breakdown of a disordered conductor under an applied current or voltage drop (often thought of as a jumping 
board to fracture problems) is yet another illustration. In a simple formulation Q , it may in fact be mapped to the 
driven disordered Ising model and the graphical scheme presented here is readily applicable to it. 

The ability to solve such basic models serves as a first step toward the study of richer, more realistic physics. A 
class of problems akin to breakdown phenomena relates to the depinning and transport of elastic manifolds through 
disordered media ]l5| , |lr| . The many variations on this paradigmatic theme include charge density waves |l6| , surface 
growth J2j|[ll| , contact lines [ p6p7| , p^ ] , polymers (ll| , flux lines in superconductors |l5|,|l6| , fractures , and earth- 
quakes |16||. Another class of related examples pertains to social dilemmas, as illustrated by financial markets p8|,p9||, 
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traffic flow, and internet congestion problems [gOj . Given this wide range of interesting systems, it is worthwhile 
to acquire an elementary understanding and intuition from some simple related models. A rough "conceptual web" 
involved in the three examples we discussed in some detail is illustrated in the following diagram. 



driving force 

\ 



coupling 
(restoring force) 
I 

breakdown phenomenon 



disorder 
(pinning force) 

/ 




V 




. T 

graphical scheme 



The main gap between our discussion and a more complete description of reality is the absence of any dynamics 
fjlf in our formulation. It takes some time for an elevator rope to break, during which stress is redistributed. At a 
boarding gate, sitting people are watching, thinking, and deciding while an impatient passenger walks up to the line. 
A driven magnet responds with some delay, and different parts of the system are continually changing while they 
influence one another. By contrast, our formulation proposes algebraic equations to describe the system rather than 
equations of motion. In the spirit of effective theories, we can extend Eqs. (|2|) and (|3|) to include time dependences 
by letting s(r) depend explicitly on time and replacing the vanishing right-hand-side by a series such as 



This modification introduces correlations in time and leads to a wealth of results relating to fluctuations in time [|3l] 
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